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Abstract 

This paper presents an approximate method for obtaining truncated 
balanced realizations of systems represented by non-rational 
transfer functions, that is infinite dimensional systems. It is based 
on an approximation to the Hankel operator. 
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Introduction 


The benefits of obtaining accurate low order models are obvious to 
the control system designer. In particular, if an accurate low order 
model can be obtained, then a low order controller can be designed 
which will hopefully maintain the designed closed loop robustness 
properties when applied to the real system. Much effort has been 
given to the model reduction problem, particularly in the finite 
dimensional situation. While input-output methods were prevalent 
in the 1960’s and 1970’s [1], the 1980’s saw tremendous effort in 
the state space approach to model reduction [ 2 ]. 

Infinite dimensional systems, on the other hand, have always been of 



interest, but have not always been so easy to deal with. Part of this 
problem comes from the extensive training that the typical control 
engineer receives in finite dimensional systems via rational 
polynomial transfer functions. The existence of a non-rational 
transfer function, other than a time delay operator, is often found to 
be surprizing. Furthermore, the origins of such a transfer function 
are usually very mysterious. This is probably due to the fact that 
control engineers typically do not understand spatial boundary value 
problems. None-the-less, non-rational transfer functions for 
describing partial differential equations have been in use since 
Heaviside and the interested reader should consult the wonderful 
book [3]. 

Many non-rational transfer functions arise from systems described 
by wave equations. These will usually have several time-delay 
operators occupying terms in both the numerator and denominator. 
Historically, even in the 1950‘s, these were being approximated by 
Pade approximations [4]. Here the individual time-delays are 
replaced by an appropriate order Pade approximation. Then after 
some tedious algebra, a rational approximation results. The 
accuracy of the resulting approximation depends significantly on the 
order of the Pade approximations used with the cost of improved 
accuracy being increased system order. 

Other than some exhausting nonlinear optimization techniques, little 
significant progress in the reduction of infinite dimensional 
systems occured until the late 1980's. This work stemmed from the 
research in the robust control area and deals with balanced 
realizations in state space and well as Hankel norms [5]. These 
methods will be discussed at length in the next section. Several 
examples will be included to demonstrate the utility of the 
approximate method developed here. 


Approximate Truncated Balanced Realizations 

Truncated balanced realizations for finite dimensional systems are 
well known and are readily implemented [2]. Many variations are 
available and the reader is encouraged to consult the large amount of 



literature for more information. Recently, Glover et al [5J have been 
able to develop this technique for infinite dimensional systems. The 
approach is similar to the finite dimensional situation in that it is 
based on singular value decomposition. In the finite dimensional 
case, the balanced realization comes from the leading terms in the 
singular value decomposition of the necessarily finite rank system 
Hankel matrix. Similarly, in the infinite dimensional case, the 
balanced realization is based on an orthogonal expansion (singular 
value decomposition) of the Hankel operator 


(Tu)(t) = 



h(t+s)u(s)ds 


where h(t) is the system impulse response. The difficulty here is 
that the Hankel operator effectively has infinite rank as nonzero 
singular values are infinite in number. The real problem in the 
single-input-single-output case is that this is an integral equation 
requiring orthogonal expansion of its symmetric kernel, h(t+s). This 
approach is referred to as that of Hilbert-Schmidt. A good 
discussion of these matters can be found in [6]. This orthogonal 
expansion for symmetric kernels is based on the eigenvalue- 
eigenfunction problem 


u n (t) = A n 


V„(S) = A n 



h(t+s) u n (s) ds 


h(t+s) v„(t) dt 


where u n (t) and v n (s) are the orthogonal eigenfunctions, or Schmidt 
pairs, corresponding to the eigenvalue A n . Once these are found, the 
Hankel kernel can then be written as 


oo oo 

h(t+s) = £ Un ^ yp( s ) = £ a n u n (t) v„(s) . 
n=l A n n=l 

The kernel can then be truncated after the least significant singular 
value cr k to yield a rank k expansion. Once this is done, it is shown 



[5] that a state space system of rank k can be developed which 
approximates the original system, with error related to (j k+ i. There 
an output normal realization is given and is repeated here; 

B; = ciiVi(O) 

C; = Ui(0) 

Ah = - 2 u f(°) 

A _ B^Bj - <TfC*Cj 

OT-OJ 

Thus, once the singular values and the Schmidt pairs are found, it is 
then a fairly simple task to obtain a truncated balanced realization. 
It should also be noted that this reduces to the standard approach 
for finite dimensional systems. Glover [5] then goes on to give the 
optimal Hankel norm approximation, but that is not pursued at this 
time. 

The problem with the approach is that it is extremely difficult in 
most cases to analytically obtain the required orthogonal expansion, 
that is, the Schmidt pairs. Thus some approximate method for doing 
this is required. The approach taken by Fredholm in the solution of 
integral equations [6] is particularly appealing for performing this 
task. Basically, the Fredholm approach is to divide the integration 
limits for the Hankel operator into evenly spaced sections in s and t, 
and then to add up the resulting sampled kernels as follows using 
the rectangular rule for integration; 


u(s) = X 



h(t+s)u(t)dt 


then for a symmetric kernel and a timestep T 


u(0) - AT[h(0+0)u(0) + h(0+ l)u(l) + • • • + h(0+n)u(n)] = 0 
u(l) - AT[h(l+0)u(0) + h(l + l)u(l) +••• + h(l+n)u(n)] = 0 

u(n) - AT[h(n+0)u(0) + h(n+l)u(l) + • • • + h(n+n)u(n)] = 0 


Clearly this is still an eigenvalue problem, but rather than finding 



Schmidt pairs which are continuous functions of time, the Schmidt 
pairs will be eigenvectors representing approximate sampled 
versions of the eigenfunctions. The above equation can be written in 
matrix form as 


1- AThoo 

-AThoi 

-XThon 


-AThio 

1-XThn 

-XThin 

Uo 

Ui 




_U n 

-ATh n o 

-XThni 

• • • 1- XTh nn _ 



or 


I-XTHjru] = 0 


where H is now the Hankel matrix formed by sampling the system 
impulse response every T seconds. It should also be remembered 
that H is symmetric. Thus finding its eigenvalues and eigenvectors 
is equivalent to a singular value decomposition. Hence, the left and 
right singular vectors of H will approximate sampled versions of the 
Schmidt pairs of the orthogonal expansion of the continuous Hankel 
kernel. Its singular values however are approximations of the 
singular values obtained from the orthogonal expansion. 

Furthermore it should be noted that the singular values of the Hankel 
matrix must be multiplied by T in order to obtain the above 
approximation. Since the relative magnitudes should not change, it 
is then necessary to divide the matrices of right and left singular 
vectors by SQRT(T). As T gets smaller and smaller, then the 
approximate version converges to the exact version. The major 
problem with the approach is in running out of computer memory to 
form H as T is allowed to get small, however reasonable results can 
be obtained for fairly large T. Some examples of this are now 
presented for comparison with the actual solutions. 

Example 1: Approximation of the unit triangle impulse response 

[5,7]. This is included for comparison with the exact solution. 

Using T=.02 yeilds a 50x50 Hankel matrix and the following second 
order system using the output normal realization given above 



H(s) = .9778s + 3.4477 

s 2 + 3.6918s + 7.0187 


as compared with the actual solution taken from [5] 

H( s ) = -95 61s + 3.6402 
s 2 + 4s + 7.6145 ' 

Comparison of the impulse responses is found in Figure 1. 


To further demonstrate the process, the solution process for T=.2 is 
now presented. The sampled Hankel matrix is 

" 1 .8 .6 .4 .21 
i .8 .6 .4 .2 0 
H = .6 .4 .2 0 0 
.4.2 0 0 0 

.2 0 0 0 0 J 

The SVD is then performed on this, followed by the T-scaling to give 


0.2a = [.4036 .0647 .0247 .0145 .0109] 

7 1.6321 1.1675 .8121 .5051 .2424 1 

i 1.2105 -.2230 -1.0630 -1.2828 -.8422 j 

-~= .8050 -1.1675 -.8121 .7707 1.3174 ! 

.4435 -1.3053 .9079 .7516 -1.3078 

L .1618 -.7215 1.3140 -1.3930 .8867 J 

-^-= [ui -u 2 u 3 -u 4 u 5 ]. 

\.2 


A second order approximation is chosen, then using the output 
normal realization given above 



.6587 
-.0783 . 


C = [1.6321 1.1675] 

A _[ -1.3319 -2.2807] 
L .3752 -.6815 J 


which is given in input-output form as 



_ 1. 1665s + . 85 15 

s 2 + 2.0134s + 1.7634 ' 

It can be seen from Figure 2 that the response is not all that bad 
considering the grossness of the approximation. 

Example 2: Approximation of the unit block impulse response [7]. 

This was done using T=.1 . The resulting system is 

_ .8605s + 3.0822 

s 2 + 1.8523s + 3.4681 

with corresponding response in Figure 3. Although it doesn't look 
much like a block, the response is almost identical to that given by 
[7J. The response from a tenth order approximation is also given and 
looks pretty good. 

Example 3: Approximation of the impulse response sin(t)/t . This 

system is perhaps even more difficult than the others since it is of 
the infinite impulse response type. Using T=.05 the approximation 
obtained was 

H( S ) = 1.0169s 2 + 1.0532s + 1.0793 

s 3 + 1.2004s 2 + 1.1813s + .6198 ' 

Although the response is not bad, it did not appear to really be that 
good, as can be seen from Figure 4. This is probably due to the fact 
that this system has an infinite impulse response. The reduction of 
such systems would appear to require a great deal of memory in 
order to obtain an accurate approximation. 

Example 4: Approximation of a simplified supersonic inlet [see 8 for 
background]. Here the system to be approximated is 


H(s) = 


50 e~- 015s 
s + 2 + 50 e- 02s ' 


The approximate truncated balanced realization is found from an 
approximate impulse response which is obtained from the following 
discrete model of this system (it uses z-inverses for the delays and 
Tustin for the s term), 



H(z) = 


50 (z 6 + z 5 ) 

2002z 21 - 1998z 20 + 50z + 50 


T=0.001 . 


The approximate truncated balanced realization is 


. _ 1000 (-.0025s 2 -.3940s + 523.7) 

x I ( S ) ***■ * 

s 3 + 105.7s 2 + 8644s + 481446 

while a Pade(1,1) approximation for the delays [9] yielded 


. 50 (-s 2 + 33.33s + 13333) 

s 3 + 185.3s 2 + 12133s + 693333 

The step responses of each of these systems is plotted in Figure 5 . 
It can be seen that the approximate truncated balanced realization 
tracks better initially but that there is a little steady state error. 
This may be due to either the balanced truncation or the 
approximation to the Hankel operator. Alternatively, for 
comparison, the Pade approximation has a much worse transient 
response, while there is no steady state error. It was also found 
that higher order truncations of the balanced realization provided 
more accurate approximations. 


Conclusion 


This paper contains an approximate method for obtaining reduced 
order models of infinite dimensional systems. It is based on an 
approximation to the Hilbert-Schmidt expansion of the Hankel 
operator. Much additional work on this approximate approach is 
necessary. In particular, determining the error bounds associated 
with this approximation is of primary importance. Furthermore, 
better approximations to the integral are being considered. 
Preliminary results indicate that the trapezoidal rule can give a 
better approximation than the rectangular rule used here. Other 
methods will also be considered. Another approach is to allow 
variable values of T for infinite duration impulse responses. Finally, 
an approach to using frequency domain information or input/output 
data should be pursued as that is all is available in some cases. 
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re 1. Truncated balanced approximation for triangular pulse and its approximation 
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0 Figure 2. Exact and crude (T=0.2) approximation of triangular pulse. 
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1 n Figure 3. 2nd and 10th order approximations to unit block. 




Figure 4. Exact and approximate sin(t)/t. 
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Figure 5. Discrete, Pade, and approximate inlet responses. 
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